Phonon Density of States in Nd(Oi_a;F2;)FeAs 
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We report measurements of the phonon density-of-states in iron oxypnictide superconductors by 
inelastic x-ray scattering. A good agreement with ab-initio calculations that do not take into account 
strong electronic correlations is found, and an unpredicted softening of phonon branches under F 
doping of these compounds is observed. Raman scattering experiments lead us to conclude that this 
softening is not related to zone center phonons, and consequently imply an important softening of 
the relevant phonon branches at finite momentum transfer Q. 
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In the last few months, the discovery of a new family 
of iron-based oxypnictides superconductorsi*^ has been 
of great interest to the community, in particular because 
of their high critical temperatures, that now reach up to 
55 K in SniQ.gFo.iFeAs^. Their study is also motivated 
by the striking similarities of these new compounds with 
high- Tc cuprates (layered structure, doping-induced su- 
perconductivity, strong electronic correlations, proximity 
of magnetic phases ...), for which there is still no consen- 
sus on the origin of superconductivity despite 20 years 
of intense research. The origin of superconductivity in 
oxypnictides is already controversial. Experimentally, 
either conventional superconducting properties (e.g. a 
BCS-like gap^) or, on the contrary, occurrence of strong 
coupling values for the superconducting parameters^ are 
found. On the other hand, on theoretical grounds, seri- 
ous arguments based on the strongly correlated character 
of these materials suggest a more unconventional super- 
conductivity mechanism^. According to recent ab-initio 
lattice dynamical calculationsii^, the Migdal-Eliashberg 
parameters A and cj„ cannot explain superconductivity 
above roughly 1 K in this class of compounds. A pure 
electron-phonon coupling as a mechanism for supercon- 
ductivity is thus likely excluded, but still, as in supercon- 
ducting cuprates where a variety of unusual behaviours of 
the lattice vibrations have been reported, studying lattice 
dynamics in oxypnictides can provide valuable informa- 
tion. To date, only polycrystalline oxypnictides samples 
are available, therefore preventing any complete phonon 
dispersion studies. In this paper, we determine exper- 
imentally the phonon density-of-states (PDOS) of the 
compound Nd(0i_2^Fa;)FeAs (x = 0.0 and 0.15) using 
inelastic x-ray scattering (IXS). 

Polycrystalhne NdOFeAs and NdOo.ssFo.isFeAs sam- 
ples were prepared by a two-step solid state reaction 
method. Stoichiometric amounts of NdAs (prepared by 
reacting Nd and As at 900 °C for 12 hours), Fe, Fe203 
and NdF3 were homogenized using mortar and pestle and 
pressed into pellets using an 8 mm steel dye and a 10 ton 
press. The pellets were inserted in tantalum cells and 
then heated in vacuum sealed quartz tubes for 2 x 24 



hours at 1150 °C with an intermediate homogenization. 
X-ray powder diffraction measurements were performed 
at A = 0.7 A on the Swiss Norwegian Beamline at the 
European Synchrotron Radiation Facility (ESRF). The 
GSAS suite of programs and EXPGUI graphical user in- 
terface were used for Rietveld fitting^ii^. The results are 
shown in Figure [1] No impurity phases were found in 
the parent compound while 5% of the NdOF phase were 
found in the doped sample. 

The IXS experiment was carried out on beamline ID28 
at the ESRF using the silicon (999) configuration at 
17.794 keV with an instrumentation energy resolution 
of 3.0 meV. Within the limit of single phonon scattering 
(Stokes process), the intensity of the inelastically scat- 
tered signal with energy (momentum) transfer uj (Q) is 
directly proportional to the dynamical structure factor—: 

j n 

jl + n{u!Qj,T))S{uj - ujQj) 

where /„(Q) is the atomic form factor of the n*'' atom 
(at position rn) in the unit cell, Af„ its mass and Wn the 
Debye- Waller factor. The summation index j refers to 
the j^^ phonon branch and o'„(Q, j) and ojqj to the asso- 
ciated eigenvector and eigenvalues, respectively. Finally, 
n(w,T) = {exp{uj /ksT) — stands for the Bose fac- 
tor. The results presented here have all been obtained at 
room temperature. In the case of a powder, the dynam- 
ical structure factor must be averaged over the sphere 
of radius Q = ||Q||, and, as previously demonstratedi^, 
averaging several measurements over a large range of Q 
{Qmin < Q < Qmax) givcs acccss to & quantity close to 
the generalized PDOS. On ID28, the scattered photons 
are analysed by a set of 8 analysers that have been cen- 
tered on two different set of angles to cover a Q range 
from 53 nm^^ to 73 nm~^ {{Q) = 63 nm^^). In this 
case, the dynamical structure factor reduces to 
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FIG. 1: (Color online) a) measured (dots) and calculated 
(solid line) synchrotron X-ray powder diffraction profiles 
(A=0.7 A) for NdOFeAs at room temperature. The lower 
solid lines show the difference profiles and the tick marks 
show the reflection positions, b) measured (dots) and calcu- 
lated (solid line) synchrotron X-ray powder diffraction profiles 
(A=0.7 A) for NdO(0.85)F(0.f5)FeAs at room temperature. 
The lower solid lines show the difference profiles and the tick 
marks show the reflection positions (from top to bottom NdO- 
(0.85)F(0.f5)FeAs, NdOF). 



u) = A{Q)^{1 + n{Lu,T)) 
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In reality, the single phonon scattering limit does not 
hold, and multiphonon contributions to the IXS response 
has to be taken into account. Although it is clear that 
the recursive procedure described in refs. [l^fisl is only 
exact in the case of monoatomic systems, it enables us 
to subtract most of this multiphonon contribution from 
our data. 

In the upper panel of Figure [2l we show the quan- 
tity lineiasuci^) X [ll> / {1 + n{uj , T))] (dashcd lines), where 
the inelastic signal hneiasUci^) is obtained after sum- 
mation of individual spectra weighted by the analyser 



efficiencies and subtraction of the elastic contribution 
(see inset of Figure [5]) , as well as the generalized PDOS 
obtained after multiphonon subtraction for pure Nd- 
OFeAs. The broad feature extending up to 65 meV in 
hneiasticiio) X [uj / (1 + n([j,T))] is essentially removed 
when removing the multiphonon contribution, and in 
the remaining spectra, we observe essentially 3 peaks of 
about 5-6 meV full-width-half-maximum located around 
13, 24 and 32 meV. In the lower panel of Figure [2l we 
report the theoretical results obtained by Boeri et al4 on 
a similar compound, LaOFeAs. To allow a direct com- 
parison between the experimental data and the theoret- 
ical calculation, we need to weight each calculated par- 
tial density-of-states G„(a') by its corresponding factor 
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/„((Q))2M-i (see tableH. 
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FIG. 2: (Color online) Comparison of the generalized phonon 
density-of-states of NdOFeAs (upper panel) with the results 
of a density functional pertubation theory calculation (lower 
panel - see details in the text). The dashed curve in the 
upper panel corresponds to the quantity Ii„eiastic{i^) x [a;/(l-|- 
n{ui,T))], before multiphonon contribution subtraction. The 
areas have been normalized to one. The inset of the upper 
panel shows the IXS measurement prior to the subtraction of 
the elastic line. 

The Debye- Waller factors we used have been deter- 
mined experimentally by high resolution x-ray diffrac- 
tion on samples from the same batch (see details in ref. 
IT^ - note that in first approximation, we have chosen to 
apply to La in the calculation the factor experimentally 
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TABLE I: Weighting factors for the partial density-of-states 



Atom 


reL /„((Q» 


M„ 


2Wr.=U,so{Q? 


Weighting factor 


La 


30.3 


138.9 


0,25798 


5.107 


Fe 


11.46 


55.8 


0.2818 


1.7756 


As 


16.95 


74.9 


0,31752 


2.792 


O 


2.335 


15.99 


0,31752 


0.248 



obtained for Nd). The result convoluted by the experi- 
mental resolution is plotted in the lower panel of figure [H 
as well as all the partial contributions. The agreement is 
fairly good since each of the three main peaks is repro- 
duced within an error of 2 meV, that can be partially at- 
tributed to the difference between the rare earth element 
of the respective study (Nd versus La) . Qualitatively, our 
data also seem to be in agreement with previous inelastic 
neutron scattering measurements^^. 

In Figure[31 we show the evolution of the PDOS of Nd- 
OFeAs on substituting O by F. The overall shape of the 
spectra is conserved, except for the second peak. Indeed, 
we notice a striking change below 24 meV, with the ap- 
pearance of an extra contribution to the PDOS around 
21 meV. In order to rule out any effect of the energy res- 
olution or the Q-sampling, we recorded the relevant part 
of the PDOS with a higher energy resolution of 1.8 meV 
(using the silicon (11-11-11) configuration at 21.747 keV 
(not shown here). The effect was still clearly visible. 
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FIG. 3: (Color online) Evolution of the PDOS of NdOFeAs 
under F:0 substitution. The area of both spectra have been 
normalized to unity. The blue line corresponds to the differ- 
ence between the two spectra. 

Simulations of the substitution of O by F using the vir- 
tual crystal approximation (VGA) show that the PDOS 
obtained for the doped sample should a priori not dis- 
play any significant differences with respect to the one 
of the parent compound^. From the calculation, three 
modes are composing the middle band extending from 
20 to 28 meV. These three modes are Raman-active, as 



reported by Hadjiev and collaborators on undoped RE- 
OFeAs (RE = (La, Sm))^^. We emphasize that recent 
measurements in CeOFeAsi^ and in NdOo.82Fo.i8FeA9i^ 
clearly show the absence of any infrared-active features 
between 110 cm~^ (13.6 meV) and 250 cm^^ (31 meV). 
We thus naturally attribute the 21 meV bump in the 
doped spectra to a softening of at least 2 meV of one of 
the three Raman-active branches on doping. Following 
Rcf. [l^, these modes are mainly be attributed to respec- 
tively out-of phase c-axis polarized vibrations of the rare 
earth (with the Ai^ symmetry). As (Aig symmetry) and 
Fe (Big symmetry). 

In order to gain further insight, and to see if one 
can explain the previous observation by the softening 
of the Raman-active zone center phonons, we have per- 
formed a (unpolarized) micro-Raman scattering exper- 
iment on various individual grains of the powders of 
doped and undoped samples. The experiment was car- 
ried out under a microscope (xlOO magnification) at- 
tached to a LABRAM Jobin-Yvon Raman spectrometer 
(1800/mm gratings) operating in backscattering geome- 
try, at 514.52 nm with an incident laser power density 
of about 5 X 10^ W/cm?. The results shown in Fig- 
ure |4] are in full agreement with those previously ob- 
tained by Hadjiev and collaborators^^ on LaOFeAs and 
SmOFeAs. We can clearly distinguish 2 peaks in the Ra- 
man response of the parent compound. The first one is 
located around 20.5 meV (166 cm~^) and corresponds to 
the Aig mode of Nd. This frequency is close to the one 
expected from the dependence of the phonon frequency 
on atomic masses and nearest neighbor distance given in 
Refs. |llli3. The distance Nd-0 is found to be 2.3149 
A in pure NdOFcAsi^. The expected renormalization 
of the rare earth Aig phonon with respect to LaOFeAs, 

given by ujNd/uJLa ^ \/imLadl^_Q/{mNdd%^_Q)) ~ 
1.009, is very close to the experimentally observed one 
^Nd/^La =1.022. The same approach for Sm gives 
^Nd/^Sm = 0.984 for an observed value of 0.976. The 
second peak around 25 meV (200 cm~^) is composed of 
two contributions^^ respectively arising from the As Aig 
mode (24.5 meV/197 cm~^) and from the Fe Big mode 
(25.3 meV/204 cm'^). As 15% of O is substituted by F 
in NdOFeAs. an increase of dNd-o of about 0.0113 A is 
founcU^. Following the above arguments, this would lead 
to a very small softening of ~ 0.15 meV of the Nd Aig 
phonon, in good agreement with the experimental finding 
of a softening of about 0.3 meV for this line. Similarly, 
a small softening of the two other branches (^ 0.3 meV) 
would be expected from the small increase of dpe-As m 
the doped compound but experimentally, a hardening of 
about 0.6 meV of these two branches is found. 

According to recent dynamical mean field theory 
(DMFT) calculations on LaOFeAs, the parent compound 
behaves like a bad semiconductor, with a small gap in the 
electronic density-of-states, separating the quasiparticle 
peak from the Fermi level^. Upon F-doping, simulated 
in the VGA, the quasiparticle peak gets closer to the 
Fermi level, and the doped compound becomes more and 
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changes in lattice parameters. On the other hand, the 
zone center softening of the Nd Aig mode seems consis- 
tent with these changes, but may be too small to be held 
responsible for the softening effect observed in the PDOS 
under F:0 substitution. The latter is certainly related 
to electron-phonon coupling in oxypnictides and might 
then be attributed to a phonon softening at finite mo- 
mentum transfer Q. Further experiments on single crys- 
tals, as well as lattice dynamics calculations taking into 
account strong electronic correlation effects, are thus re- 
quired to understand this softening in the PDOS. Com- 
bined together, they will shed new light on the impact of 
electron-phonon coupling on the already rich physics of 
oxypnictides. 



FIG. 4: (Color online) Raman scattering spectra of NdOFeAs 
(black points) and NdOo.ssFo.isFeAs (red points). The dot- 
ted lines are the individual fits for each of the three phonon 
lines and the solid line is the sum of the fitted lines. 

more metallic. As a consequence. Coulomb screening is 
modified and one can expect variations of the frequency 
of phonon branches coupled to the electronic continuum. 
This effect will be likely most pronounced for the phonon 
involving the motion of atoms of the conducting FeAs 
layer since most of the states close to the Fermi level 
originate from the Fe 3d bands. This could possibly ex- 
plain the zone center hardening of the FeAs layer phonons 
found in our Raman measurement, unexpected from sole 



Note 

After first submitting this paper, we became aware of 
similar PDOS ncutrorii^, x-ray^S^ as well as nuclear rcs- 
onani*2^ inelastic scattering measurements on iron based 
oxypnictide superconductors. 
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